Effect of electron interactions on the conductivity and exchange coupling energy of 

disordered metallic magnetic multilayer 
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We consider the effect of electron-electron interactions on the current-in-plane (CIP) conductivity 
and exchange coupling energy of a disordered metallic magnetic multilayer. We analyze its de- 
pendence on the value of ferromagnetic splitting of conducting electrons and ferromagnetic layers 
relative magnetizations orientation. We show that contribution to the CIP conductivity and ex- 
change coupling energy as a periodic function of the angle of magnetizations relative orientation 
experience 2n — > n transition depending on the characteristic energies: ferromagnetic splitting of 
the conducting electrons and the Thouless energy of paramagnetic layer. 
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I. INTRODUCTION 

The most interesting features of perfect metallic mag- 
netic multilayer are the oscillatory behavior of bilinear 
exchange coupling energy between ferromagnets^ 2 - due to 
Friedcl oscillations, and large magnetoresistance in small 
magnetic fields 3 -! 4 -. 

The magnetic structure of adjoining magnetic lay- 
ers in perfect multilaycred structure oscillates between 
ferro and antiferro-magnetic states with increasing the 
spacer thickness L. Disorder in layers contributes to bi- 
quadratic exchange coupling. It was shown that fluctu- 
ations of thickness of the paramagnetic layer give rise 
to biquadratic exchange coupling, often leading to non- 
collincar magnetic ordering^. In Ref. 00 the role of scat- 
tering of conducting electrons by impurities in metallic 
magnetic multilayer was studied. It was pointed out that 
in the case of small (compared to thickness of the layers) 
mean free path of conducting electrons, when average 
Friedel oscillations are exponentially suppressed, the ex- 
change coupling energy due to random Friedel oscillations 
and correlation effects can have biquadratic form. Tran- 
sition to the non-collinear phase in disordered structure 
with increasing L was experimentally observed in Ref. 8. 

It is established that the magnetoresistance of per- 
fect metallic magnetic multilayered structure is related 
to spin-depended scattering of conducting electrons at 
the interfaces between the layer o^ 10 ' 11 . In the case of 
disorder, scattering of conducting electrons by impuri- 
ties suppresses the effect of spin-dependent scattering on 
magnetoresistance. Moreover, it was theoretically shown 
that the magnitude of the magnetoresistance decreases 
exponentially in the case of the current-in-plane (CIP) 
geometry, when thickness of the nonmagnetic spacer ex- 
ceeds the mean free path of conducting electron a 10 i 12 . 

In the present paper we study the effect of electron- 
electron interactions on CIP conductivity and exchange 
coupling energy of disordered metallic multilayered struc- 
ture consisting of two ferromagnetic layers with param- 
agnetic spacer. We consider the case when mean free 
path of conducting electrons is smaller than thicknesses 



of layers. It is known that in the disordered conductors 
electron-electron interactions result in anomalous contri- 
butions to conductivity, thermodynamic quantities and 
negative magnetoconductivit y 13 > 14 . Physics behind the 
effect of electron-electron interactions in disordered con- 
ductors is the electron scattering by the random Friedel 
oscillations 1 ^. Frcidel oscillations in magnetic multilay- 
ered structure do depend on the angle tp between direc- 
tions of magnetization in magnetic layers. The study of 
this dependence is the subject of this paper. 

The scattering by Freidcl oscillations exists in any type 
of ferromagnetic structure: itinerant (d- type) or local- 
ized (/- type). In what follows, the most relevant fac- 
tors of disordered magnetic multilayered structure are 
the characteristic of the disorder - the Thouless energy 
D/L 2 (D is the conducting electron diffusion constant), 
and ferromagnetic splitting of conducting electrons. Be- 
cause of that we consider the model of localized s-f mag- 
netism that contains these parameters in the most trans- 
parent way. We propose that each ferromagnetic layer is 
described by the homogeneous magnetization, and that 
directions of magnetizations of different layers make an 
angle ip. 

The conductivity and exchange coupling energy in 
magnetic multilayered structure are periodic functions of 
(p. We show that depending on the ratio of introduced 
above characteristic energies, the angle ip dependent con- 
tributions of electron-electron interactions experience the 
2tt — > 7r periodicity transition. Magnitudes of contribu- 
tions are estimated. 

We suppose that our results might be relevant to the 
series of works related to transport properties of meso- 
scopic fcrromagncts with domain walls, where it was 
cxperimentalljii^ and theoreticall y 17 ' 18 shown that effects 
of electron interactions and weak localization are impor- 
tant. 

We would like to mention that results of this pa- 
per complement the theory of electron-electron interac- 
tions in disorders conductor o 13 ' 14 as a study of effect of 
electron-electron interactions in spatially inhomogeneous 
effective magnetic field that is imposed by the magneti- 
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zation in layers. 



II. THE MODEL 

The system of our study is disordered metallic mag- 
netic multilayer structure. It consists of two ferromag- 
netic layers of the /-type located at L/2 < \z\ < L\/2 and 
a paramagnetic layer between (\z\ < L/2). The Hamil- 
tonian of conducting electrons of this system is: 

H = H + H c +H sf . (1) 

Here H is the Hamiltonian of free electrons in a ran- 
dom field V (r) . Treating scattering of free electrons we 
use the standard diagram technique 155 ., which assumes 
(V(r)) = and (V(r)V(r')) = 1/(2-kv t)5 (r - r'). 
Here vq is the density of states at the Fermi level per 
one spin, r is the electron mean free time, and we have 
set h = 1. We assume the Hamiltonian of free electrons 
to be the same in every layer. 

The second term He is the Hamiltonian of the 
Coulomb interaction between conducting electrons. 
H c = ±J*i(r)^(rO^*Mr0^r)drdr', where 

^J, (r) and ^ (r) are the electron creation and annihi- 
lation operators correspondingly. We will treat Coulomb 
interaction within the random phase approximation. 

The third term H s f is the Hamiltonian of s — f ex- 
change in the ferromagnetic layers. At temperatures 
much lower than the Curie temperature one might ne- 
glect electron-magnon interaction, therefore 

i 

ISn s JdT^ a (r)(n (r)-a^(r) (2) 

F 

Here Si is a spin of localized /-electrons, ns is their den- 
sity. / is the s — f exchange interaction, n (r) is the 
ferromagnetic layer magnetization direction unit vector. 
Integration here is over the ferromagnetic layers. 

Neglecting the contribution of electron-magnon inter- 
action to the conductivity and exchange coupling energy, 
we assume that the Coulomb energy per electron is larger 
than the ferromagnetic splitting ISng. 

Finally we consider the following Hamiltonian of a dis- 
ordered metallic magnetic multilayer: 

H = H +H C + ^ J dr *t (r) (n (z) . fj a ^ (J (r) . 

F 

(3) 

This Hamiltonian describes a Fermi liquid in an effective 
inhomogeneous magnetic field that acts only on electron 
spin. This magnetic field results in the spin splitting with 
energy e exc = 2ISns in the ferromagnetic layers, and in 
zero splitting in the paramagnetic layer. Direction of 
this magnetic field is defined by the ferromagnetic layer 
magnetization direction unit vector n(z): n (z) = ni at 
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FIG. 1: A. Feynman diagram representation of the Hartree- 
type corrections to conductivity. The solid lines in the dia- 
grams denote the disorder averaged Green's functions, wavy 
line presents the screened Coulomb interaction, dotted line 
with the cross represents impurity scattering, and the black 
square stands for the diffusion ladder. B. Hartree-type cor- 
rection to thermodynamic potential. C. Equation for the dif- 
fusion ladder 

—Li/2 < z < -L/2 and n(z) = n 2 at Li/2 > z > L/2, 
with (ni • 112) = cosip. 

Now we are ready to study the effect of electron inter- 
action on the CIP conductivity and exchange coupling 
energy between the ferromagnetic layers. We will be in- 
terested in its dependence on the angle between the direc- 
tion of magnetization tp, ferromagnetic splitting energy 
e exc , and paramagnetic thickness L. 

A. Conductivity 

We use Kubo linear response formalism in order to 
study the effect of electron interaction on conductivity. 
Treating electron interaction as a perturbation, we come 
up with corrections to conventional Drude formula for 
conductivity. We study only the first order perturba- 
tion correction. The magnetic field dependent, or in our 
case magnetization dependent, correction to conductivity 
comes from the interaction between two electrons in the 
triplet state (so-called the Hartree-type correction) 13 ' 14 . 
Disorder averaged correction can be described with the 
help of a Feynman diagram, shown in the fig. 1A. Black 
squares in the figure stand for the diffusion ladder that 
represents electron's charge and spin densities propoga- 
tion. Calculations of presented diagrams for the homoge- 
neous magnetic field are carried out in Ref . [l3|. We make 
use of the bilinear representation of the diffusion ladder 
in order to calculate these diagrams for inhomogeneous 
magnetic field. Our calculations show that in the case of 
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the CIP geometry all diffusion ladder combinations that 
appear in the figure can be simplified in to the following 
expression: 

6a(z)=iF(z)^ 

OO 

x / dw^(wcoth(^,))Trl^°(at,at J 0,-iw) ) (4) 

— OO 

where F (z) = F fer 6 (\z\ - L/2) + F par 6 [L/2 - \z\) is 
the electron interaction constant that is discussed in 
Ref. [ID. We assume that interaction constants might 
be different in ferromagnetic and paramagnetic layers. T 
is the temperature, T>"P is the diffusion ladder. For the 
CIP geometry it is convenient to consider conductance 
SG = J dzSa (z) as a sum over the paramagnetic and 
ferromagnetic layers. 



B. Exchange coupling energy 

Contribution to the multilayer thermodynamic poten- 
tial, which depends on the relative orientation of the fer- 
romagnetic magnetizations, is also due to interaction be- 
tween two electrons in the triplet state. This correction 
is presented in the fig. IB and given by an expression^ 



<$n(y) _ 

s ~ 



(5) 



T £ \u n \f4$ r fdzF(z)Tr>DZZ(z,z ) q,u, n ), 

w„|t<1 



where u>„ = 2~imT is the Matsubara frequency, S is the 
area of the system. The angle-dependent part of 8il (tp) 
determines exchange coupling energy of adjoint ferro- 
magnetic layers. 



C. Diffusion ladder 

The angle tp dependence in expressions ([4]) and (O is 
due to the diffusion ladder. Graphical equation for the 
diffusion ladder T>"P (z, z', q, ui n ) is shown in the fig. 1C, 
and in the case of diffusion approximation (|a; n |T < I, 
£exc T < 1) and \q\vpT < I with vp being the Fermi ve- 
locity) satisfies the following differential equation 



pa/3 



(-D£ + Dq2 + \u 
n (z) (cr Q7 X>^ - V^(T 1V ) signw n = (6) 



5(z - z')d a p5,j, v , 



Here D = \v f t is the diffusion constant. We assume 
it to be the same in every layer. We notice that in 
expressions (j4|) and ((5|) only trace of the diffusion lad- 
der appears. Since the trace is invariant under unitary 
transformations, we may choose the ferromagnetic lay- 
ers magnetization direction unit vector in the convenient 
for further calculations: in — L\/2 < z < —L/2 the di- 
rection is rix = (1,0,0), and in Li/2 > z > L/2 it is 
ri2 = (cos ip, 0, sin tp). 



Solution of the equation ([6]) has to satisfy the boundary 
conditions -j^L)^ (z, z' , q, ui n ) = at z = ±Li/2, and the 
continuity conditions at z = ±L/2. 



III. 



2tt 



tt TRANSITION 



Here we present our main results by skipping all deriva- 
tions that will be given in the appendix to the paper. We 
consider only angle-dependent parts of j4| and ([5|). 

Characteristic energies of a magnetic multilayer are the 
Thouless energy of the paramagnetic layer D/L 2 and fer- 
romagnetic splitting energy e exc . We present results for 
the cases of large ferromagnetic thickness d — L\/2—L/2: 
e exc d 2 /D > 1 and for small temperatures: e exc >> T. 

Our main finding is the following: depending on the 
ratio of the Thouless energy and ferromagnetic splitting 
energy, contribution to the CIP conductance ([!]) and ex- 
change coupling energy between ferromagnetic layers ([5]) 
might have 27r or 7r periodicity as a function of ip. 2n — > tt 
transition occurs when ratio e exc L 2 /D increases. 

(i) In the case of large Thouless energy e exc L 2 /D < 1 
the leading term the in angle-dependent parts of the CIP 
conductance and exchange coupling energy is propor- 
tional to cos ip. One can show that ferromagnetic layers 
give the main contribution to the conductance as well as 
to exchange coupling energy: 



sn(ip) 



5G 



, COS ip 



-F 



fer 



32tt 2 



COS tp, 



(7) 



2(4tt)^ D 



0.06F par + 0.7F fer In 



D 



,L 2 



(8) 

Correction to conductance originating from the param- 
agnetic layer is reduced by a factor of */ e exc L 2 / D . 

(ii) In the limit of e exc L 2 /D > 1 the leading term in 
(HI) and © is cos 2 tp, and it is independent on the ferro- 
magnetic splitting energy e exc : 



SG~- (0.7F fer + 0MF par ) cos 2 tp, 



Sft (tp) (3F fe 



s 



2(4tt) z D 



cos 2 tp. 



(9) 



(10) 



The physics behind this transition is the following. Ac- 
cording to Ref. [3 the magnetic field dependent contribu- 
tions to the conductivity and exchange coupling energy 
are due to interaction of electrons with total spin J = 1 
and projection J z = ±1. In our case when direction of 
n (r) varies in space, projections J z = 0, ±1 are defined 
only locally. 

At small e exc the leading contribution is due to inter- 
action of electrons with total spin J = 1 and projections 
J z = ±1, and the dependence of the conductivity and 
of the exchange coupling energy is proportional to cos tp. 
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In the case of large e exc diffusion modes with total spin 
J = 1 and projection J z = ±1 do not penetrate into fer- 
romagnetic layers and their contributions to the conduc- 
tivity and the exchange coupling energy are suppressed. 
The main contribution in this case is due to the diffusion 
mode with J = 1 , J z = and proportional to cos 2 (p. Let 
us note, that in this case the electron-magnon interaction, 
which creates J z = ±1 diffusion modes in ferromagnetic 
layers might be neglected even if the Coulomb energy per 
electron is smaller than the ferromagnetic splitting. 

At temperatures T > D/L 2 the contribution to the 
conductance decreases as 

e 2 2L 2 

5G = -^2-^rl F par+5F fer } C OS 2 <p. (11) 

And exchange coupling energy decreases with the tem- 
perature as exp(—Ly / %-KT/D). 

In the Fermi liquid theory the interaction constants are 
considered as some parameters. In the case of weakly non 
ideal 3D Fermi gas F = k 2 / (2p 2 F )ln(l + App/n 2 ) , where 
k and pp are the inverse screening length and the Fermi 
momentum, correspondingly^. For the Coulomb inter- 
action F par and Ff er are both positive. In both limits, 
(i) and (ii), contribution to conductance has a minimum 
at ip = 0. Exchange coupling energy has a minimum at 
ip = in the first limit, and at (p = tt/2 in the second 
limit. 

Let us estimate the value of exchange coupling energy. 
Taking D = 10cm 2 / sec, e exc = 300A' x ks (where ks is 
the Boltzmann constant), we obtain that crossover occurs 
at the paramagnetic layer thickness of L = 'TiD / 'e exc = 
5nm. Plugging these values into equations and (ITU)) , 
and assuming Ff er ~ F par ~ I we obtain SCI (ip) / S ~ 
0.5 x 10~ 3 erg/cm 2 . 

Correction to conductance is of the order of percents 



IV. CONCLUSION 

In the present paper we have considered properties of 
a disordered magnetic metallic multilayer in the case of 
small electron mean free path compared to paramagnetic 
layer thickness. The angle <p dependence of the CIP con- 
ductance and of the exchange coupling energy is deter- 
mined by the electron-electron interactions. Depending 
on the ratio of the Thouless and ferromagnetic splitting 
energies these quantities might be 2tt as well as n periodic 
functions of ip. 

Obtaining our results we have neglected the effect of 
spin-orbit scattering. This scattering smears out the 
magnetic dependence of electron interaction corrections 
to thermodynamic and kinetic properties^. However, 
this effect become crucial for paramagnetic thickness 
(larger than the spin relaxation length due to spin-orbital 
scattering) when values of calculated quantities already 
become too small. 



We would like to point out that in Ref. contribution 
to exchange coupling energy due to the interaction was 
considered only in the paramagnetic layer. In this paper 
we have shown that the contribution of ferromagnetic 
layers is larger than that of paramagnetic layer. 
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APPENDIX A: DERIVATIONS 

Here we describe the main steps of derivation. We solve 
the differential equation ^ for diffusion ladder. Accord- 
ing to the expressions ((4]) , ([5]) the solutions of interest are 
those with z and z' lying in the same layer. 

(i) Let us first derive the diffusion ladder when z' < 
—L/2 is in ferromagnetic layer. It is convenient to present 
solution in paramagnetic layer (\z\ < L/2) in the follow- 
ing form: 

Vf» = A >~ QZ + UU B lsU Sv eQ*. (Al) 

And in ferromagnetic layer at —L\/2 < z' < z < —L/2 
solution that satisfies the boundary condition at z = 
—L\/2 is convenient to present in the following form: 

F$ = + E ± PFCr&P£ coshQ (z + if) 

+p *-y E 70 p \u cosh (q x ( z +lf)) (A2) 

+P ^M^P^ cosh (Qt {z+%)), 

where 

•^0(11/ = (A3) 
E± % P ± cxp (-Q(% + z)) coshQ + kf) 
+ P+ dq; exp (-Qi(^ + z)) coshQ! {zf + ^) + 

+ ^Sf CX P (-Qi(tt + z )) cosh Ql i z ' + ¥) . 

where P± = (1 ± o z )/2 are projectors on the spin up and 
down states correspondingly; Q = i/q 2 + |a;„| /D, Q\ = 

\/q 2 + (\u> n \ + if-exesignun) / D , and U = exp(iiptj y /2) is 
the matrix of rotation along the y axis. A, B, C, E, M 
are matrixes-constants to be determined with the help of 
continuity conditions. 

It is convenient to obtain closed set equations for ma- 
trixes A ± B with the help of the continuity conditions 
for rotated diffusion ladder JM]) V$ = U^V^U^ at 
z = L/2: 

■R(Q)P^Vf x P^ = niQ^p^vilPh" = 

UiQDPZ'VllPl" = 0, (A4) 
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where we have used the differential operator H(Q) = with£(Q) = (Qsinh( Q( ^ +Ll) ) - cosh ( 
/ „ . , / ni t.- r., > \ . f ntT.-T.,\\ j \ . V V / V 



Q(-L+Li 



(q sinh ( 



Q(L-Li) 
2 



cosh 



/ Q(L-L 1 ) \ 
^ 2 J dzj 



d z 



ity conditions for diffusion ladder (|A1|) at 



And continu- 
al = —L/2 are 



£(Q)PFV™P£ 



t{Q 1 )pr i vi{ph v 



■X0 jjfJtV 



coshQ (2'+ -^l) , 



D 

t{QDP^Vf x P^ =^tl C oshQt(z' + ^) 



coshQi (z' 



^1 
2 



Solving for matrixes A, B and subtracting angle inde- 
pendent part we obtain the trace of the diffusion ladder 
for ferromagnetic layer (|A2|) 



-4Ai 

Q 



J Tr^(z,z,q,\u n \)dz = 

\z\>L/2 

Rc(1+A !T +C ° syl + (A ^ — A, Ai -» -Ai) 



sinh(2Qrf) _|_ 



-^Re 



2+A+AA 1 +(2A 1 +AA 1 +A)cos ¥ > a a 

2Qi{Q2 _ Ql)w+ + (A -> -A, Ai 

I 



-Ax] 



(A5) 



where we have used next notations: 

a _ Qi sinh(dQi)~Qcosh(riQi) (_riT\ 
n _ Qi sinh(dQi)+Q cosh(dQi) LX P ' 

W± = Re (1 ± A*) [1 + AAi ± (Ai + A) cos tp] , 
Ai = cxp (-QLi) . 

(ii) The same procedure can be applied for the case 
when z' is in paramagnetic layer. Now the solution of 
equation ([6]) when z is in paramagnetic layer is convenient 
to present in the following form: 



In this case the continuity conditions for the rotated dif- 
fusion ladder {56) = U^VJ^U^ at z = L/2 are 
the same as (|A4p . At z = —L/2 continuity conditions for 
the diffusion ladder (|A6j) change as 



V 



C{Q\)P^VJ X P^ 



-Q\z-z'\ 

2QD 



(A6) 



Solving for matrixes A, B and subtracting angle indepen- 
dent part we obtain the diffusion ladder for paramagnetic 
layer (|A6I) , whose trace is 



/ 

\z\<L/2 



TrV^(z,z,q, \tu n \)dz 



LR.cA*(X+A) 
2DQW+ 



(COS if + A X ) + (1 + A 1 COS If) 



-AiL 



Rc[(l+A")(l+A co S y)] 5i °^ I - -Rc[(l+A'')(A+cosy)] 
2DQW+ 



LRcA'(l-A) 
2DQW- 



cos <p + Ai) + (1 - At cos ip) 

Re[(X-A*) (X-A cos tp)] aln gf L -Re[(X-A) (A* -cos tp)] 



2DQW- 



r 



(A7) 



The results (|T) — (fTTj) will be recovered if one will plug 
diffusion ladder traces (|A5j) and (|A7j) into the expres- 
sions for the conductivity @ and exchange coupling en- 
ergy ([5]). The integrals over the frequency in resulting 
expressions can be evaluated as following: in the case of 
e exc L 2 /D, e exc d /D > 1 the main contribution is due to 



frequencies \u> n \ < e exc where A ~ Re A ~ exp (— QL). 
In this case expressions (j4]) and ([5]) are even functions of 
cosip. In region from e exc L 2 /D < 1 the main contribu- 
tion is due to \u> n \ > e exc in this case Re A ~ |A| 2 , and 
expressions ((4]) and ([5]) have the 2tt periodicity. 
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